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A GENERIC MODEL FOR SPOUSE’S PENSIONS WITH A 
VIEW TOWARDS THE CALCULATION OF LIABILITIES 

ALEXANDER SOKOL 


Abstract. We introduce a generic model for spouse’s pensions. The generic 
model allows for the modeling of various types of spouse’s pensions with pay¬ 
ments commencing at the death of the insured. We derive abstract formulas 
for cashflows and liabilities corresponding to common types of spouse’s pen¬ 
sions. We show how the standard formulas from the Danish G82 concession 
can be obtained as a special case of our generic model. We also derive ex¬ 
pressions for liabilities for spouse’s pensions in models more advanced than 
found in the G82 concession. The generic nature of our model and results 
furthermore enable the calculation of cashflows and liabilities using simple 
estimates of marital behaviour among a population. 


1. Introduction 


The motivation for this paper is the accurate calculation of the liabilities corre¬ 
sponding to the particular type of life insurance policies known as spouse’s pen¬ 
sions. In such a policy, payments are made to the spouse upon the death of the 
insured, in the case where a spouse is present. Many pension funds offer products 
such as this and have a considerable interest in efficient practical estimation of 
their corresponding liabilities to the policyholders. 

In the Danish G82 concession, specifying many standard types of life insurance 
policies, several types of spouse’s pensions are described, see e.g. [3] for more 
on this. The G82 concession describes specific formulas for the calculation of the 
liability, meaning the expectation of the discounted value of the future payments, 
under certain assumptions. The formulas, however, are built around a very par¬ 
ticular and unwieldy model, for which recent estimates of the parameters are not 
generally known. Our main objective in this paper is to develop a flexible modeling 
framework for estimation of liabilities for spouse’s pensions. 


Date: March 9, 2015. 

2010 Mathematics Subject Classification. Primary 62P05; Secondary 60G55. 

Key words and phrases. Life insurance, Liability, Gashflow, Marked point process, Longevity. 


1 



2 


A. SOKOL 


Apart from this, the forthcoming Solvency II rules from the European Union has 
led to increased theoretical and practical interest in the calculation of not only the 
liability, but also the cashflow of insurance policies, meaning the expected rate of 
payments on the insurance policy in the future, see e.g. [2]. One classical setup 
for such calculations is to let e.g. the health state of the insured be modeled by 
a continuous-time Markov chain or semi-Markov chain, see For spouse’s 

pensions, the presence of a future spouse with a priori unknown age excludes the 
possibility of a simple Markov chain model, and therefore different methods must 
be applied to obtain expressions for the cashflow of such policies. Consistently with 
the above, we aim to obtain a modeling framework which enables the calculation 
of both liabilities and cashflows. 

Finally, the Solvency II rules also specify the necessity of modeling the longevity 
risk inherent in life insurance products, meaning the modeling of longevity im¬ 
provements in populations over time. In the liability expressions of the G82 con¬ 
cession, longevity improvements are not present. It is therefore of interest to 
obtain models for the calculation of the liability of spouse’s pensions in which 
longevity improvements are included, such that e.g. the mortality benchmark in¬ 
tensities with longevity improvements reported in [1] can be used when calculating 
cashflows and liabilities. 

In this article, we develop a generic model for spouse’s pensions, and derive ex¬ 
pressions for cashflows and liabilities for a wide family of pension products. We 
also show how to obtain the classical G82 concession expressions for the liabil¬ 
ity corresponding to a spouse’s pension as a special case, using a marked point 
process model. Finally, we show how to extend this model to include longevity 
improvements. 

The remainder of the article is structured as follows. In Section [51 we review 
the notions of payment processes, expected cumulative payments, cashflows and 
liabilities in the context of a simple Markov chain model. In SectionjSl we introduce 
our generic model for spouse’s pensions and derive expressions for cashflows and 
liabilities. In Section HI we show how the expressions from the G82 concession can 
be replicated in our framework through a marked point process model. Here, we 
also show how to extend this model to include longevity effects. Finally, in Section 
111 we discuss our results. Appendix fAl contains proofs. 


2. Review of the continous-time Markov chain framework 


In order to motivate our model, we first recall the modeling framework based on 
continuous-time Markov chains as discussed in e.g. [2]. Gonsider a simple life 
insurance product paying one amount of monies per time from a given timepoint 
c and onwards, for as long as the insured is alive. Let Z be the health state of the 
insured, taking the values a (alive) and d (dead). In order to model this insurance 
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product, we may consider the process 

(2.1) Bt = f l(Zt=a)l(t>c) ds- 

Jo 

For any t > 0, Bt describes the cumulative payments paid out to the insured. We 
refer to B as the cumulative payments process, or simply as the payment process. 
We may then consider 

(2.2) At = EBt = f l(^t>c)PiZt = a) ds, 

Jo 

the expected cumulative payments. Since A is continous and differentiable almost 
everywhere, we may let a denote the Radon-Nikodym derivative with respect to 
the Lebesgue measure, yielding 

(2.3) at = \(t>c)P{Zt = a). 

We refer to a as the cashflow corresponding to the insurance policy. Finally, 
introducing an interest rate model based on a deterministic short rate r, we may 
define 

poo 

(2.4) L = E 

Jo 

the liability corresponding to the insurance policy. These concepts of cumulative 
payment processes, expected cumulative payments, cashflows and liabilities, are 
well known in various guises from the literature, see e.g. [3 El 0 E]. In the next 
section, we use the same framework in the context of a generic model for spouse’s 
pensions. 


3. A generic model for spouse’s pensions 


In this section, we introduce our generic model for spouse’s pensions. We are 
interested in modeling spouse’s pensions of the type where the spouse is entitled 
to certain payments contingent upon the death of the insured as well contingent 
upon a generic ’’policy state”. Usually, this latter ’’policy state” will be the health 
state of the spouse, such that e.g. payments only are made for as long as the spouse 
is alive, but for flexibility, we do not limit ourselves as regards the nature of this 
policy state space. After the introduction of the modeling framework, we derive 
expressions for cashflows and liabilities in the generic model. Also, we illustrate 
the usefulness of our model by deriving expressions for cashflows and liabilities for 
several types of spouse’s pensions. 

Assume given a probability space Let T be a random variable taking 

its values in R+, describing the time of death for the insured. Let A be a random 
variable taking the two values Xs and Xm, corresponding to ’’single” and ’’married”, 
respectively, describing the marital state of the insured at the time of death, T. 
Let F be a random variable denoting the age of the spouse at the time of death 
T. Here, we let d ^ M+ be a ’’coffin state” held by Y if the insured was unmarried 
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at the time of death, meaning that we assume Y = d whenever X = Xs, and 
otherwise Y takes its values in K+. 

Finally, for each m, y > 0, let (Z^’'^)r>u-y denote a stochastic process on [u — y, oo) 
with some common finite state space E. We think of as a stochastic generic 
’’policy state” for the case where a spouse exists at the time of death of the insured, 
and that spouse has age y at time u. Consistently with this, we let be defined 
on [u — y, oo), where u — y is the (possibly negative) timepoint when the spouse 
then had age zero. In the most common case, Z'^'^ will describe the health of the 
spouse. Furthermore, let Dr{E) denote the space of cadlag functions from [r, oo) 
to E, and assume that takes its values in Du-y{E). Let FV denote the 
space of mappings of finite variation from ]R_|_ to M. For u,y G K+, assume given 
measurable mappings n„ j, : Du-y{E) FV, where both spaces are endowed with 
the tJ-algebras induced by the coordinate projections. We refer to the mappings 
(n„_y) as the payment functions, and define a process C by 


Ct = nT.Y.t(^^’^), 


(3.1) 


where n„ t(z) = n„ ,^(z)t. The interpretation of this is as follows. The expression 
^u,y,tiz) represents the cumulative payments made to the spouse at time t, given 
that the death of the insured occurred at time u, that the insured was married at 
that time, and that the age of the spouse at that time was y, and given the policy 
state history z since the birth of the spouse. As a consequence, C represents the 
unconditional cumulative payments for the insurance policy with payment func¬ 
tions (n„^y), excepting that C does not prescribe payments to begin conditionally 
upon the death of the insured while married. 

It remains to define the actual cumulative payment process corresponding to the 
components of the spouse pension described above, similarly to how we in Section 
[2] defined cumulative payment processes for a simple policy. To this end, we define 
a process B by 



(3.2) 


The process B then has paths of finite variation, and corresponds to the cumulative 
payment process for the spouse’s pension with payment functions (n„_y). Given 
the joint distribution of {T, X,Y, Z’^’'^), the model for spouse’s pensions is fully 
specified, and given the payment functions (n„_j,), a particular spouse’s insurance 
policy is fully specified. 

Example 3.1. We will show how the above framework contains a classical spouse’s 
pension, corresponding to the G82 concession reward no. 810. Let E = {a,d}, 
where a denotes alive and d denotes dead, corresponding to that Z^’^ denotes the 
health state of the spouse at time r > u — y, for the case where the spouse has age 
y at time u. Define, for t > 0, 


(3.3) 
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This choice of If corresponds to the spouse receiving monies at a rate of one unit 
of money per unit of time, whenever the spouse is alive, starting at the death of 
the insured at time u. To see this, note that 

(3.4) Ct = = 1(^/0) l(t>T) J l(zT’"=a) 

corresponding to Ct being equal to monies accumulated at a rate of one unit per 
unit of time, over the time period [T, t] from the time of death of the insured to 
the current time t, with payments only accumulating at time r € [T, t] when the 
spouse is alive. The presence of in the integral corresponds to the age of the 
spouse at time T being Y. Note that we need to include the indicator ioi y ^ d 
in the definition in order to ensure that the expression for If is well-defined even 
when y = d. In practice, we only take interest in the values of Ct when Y ^ d, so 
the actual values of n when y = d is not of any consequence. We can see this by 
considering the cumulative payment process of the insurance policy, which is 

(3.5) Bt= l(s>T)l(A:=a:„) dCg = l(t>7')l(x=a;„) 

Jo 

since Y ^ d whenver X = Xm- o 



As in Section!^ we may now take an interest in the cashflow and liability corre¬ 
sponding to the cumulative payment process (13.211 . In order to prove expressions 
for these, we require some regularity conditions on the joint distribution of the 
variables T, X and Y and the processes as well as regularity conditions on 
the payment mappings. 

Assumption 3.2. We assume that the distribution of T has a density h with 
respect to the Lebesgue measure on K.+ , that the conditional distribution of Y 
given T = u and X = Xm has a density f{-\u) with respect to the Lebesgue measure 
on 1R+, and that each policy state process is independent of (T, X, Y). Finally, 
we assume that for it > 0, y > 0 and z G Du-y{E), it holds that \lu,y,u-{z) = 0. 

Here, n„_y^„_(z) is the limit of Yiu^y^r{z) as r tends to u from below. The latter 
assumption that n„_y_„_(z) = 0 corresponds to assuming that no payments are 
made prior to the death of the insured. Removing this assumption would not 
complicate matters considerably, but would merely result in different factors in 
e.g. the formula ()3.6p . However, as the condition Wu,y,u-{z) = 0 holds for all 
spouse’s pensions of interest to us, we make the assumption in order to simplify 
our results. 

In the remainder of this section, we assume that Assumption 13.2| is in force. Fur¬ 
thermore, we let g{u) = P{X = Xm\T = u) for u > 0. The following three results 
yield expressions for the expected cumulative payments, the cashflow and the lia¬ 
bility of a spouse’s pension, respectively, corresponding to (12.21) , p.3p and (j2.4ll of 
the example in Section [51 
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Theorem 3.3. It holds that 

pt pCO 

(3.6) At= h(u)g[u) / f{y\u)E\lu,y,t{Z'^'^)<dydu. 

Jo Jo 


The fomula (13.61) may be interpreted as follows. When calculating the expected 
cumulative payments of the spouse’s pension, we first condition on the time of 
death while married, which has density u i—>■ h(u)g{u). Note that as this event 
does not always occur, the density u i-A h{u)g{u) is generally defective, meaning 
that it does not have unit Lebesgue integral, instead its integral over R-|_ is equal 
to P{X = Xra)- Having conditioned on the time of the death of the insured while 
married, (I3.6p yields that the expected value of the cumulative payments up to 
time t is 

^OO 

(3.7) / f{y\u)EIly,^4Z^^y)dydu. 

Jo 

This corresponds to first conditioning upon the age of the spouse at the time of 
death of the insured, which has density fi-\u) when the time of death is u. Given 
this age, the formula (13.71) then shows that the expected value of the cumulative 
payments up to time t is 

(3.8) EUy^uAZ'^’^) 

which corresponds to expected value of the payments Ily^u,t{Z'^’y) accumulated 
from time u, when the death of the insured occurred, to time t. Here, the inde¬ 
pendence of Z'^’y and (T, X, Y) ensures that the expectation does not need to be 
conditioned upon the values of T, X and Y. 

Corollary 3.4. Assume that Ay^u,u{z) = 0 and that for all y,u>0, the mapping 
t i-A EIly u,t{Z) is differentiable for t > u with a derivative which is bounded on 
compact sets. It then holds that the cashflow (at) corresponding to the cumulative 
payment process exists, and it is given by 

(3.9) fiy\u)(^-^EIlu,yAZ"'’^)jdydu. 

Corollary 3.5. Assume the regularity conditions of Corollary \3.4\ Given a deter¬ 
ministic short rate r, the liability corresponding to the payment functions {Jdy^uA 
is given by 

(3.10) L = h{u)g{u) f{y\u) dtdydu. 


In the corollaries, (13.91) is essentially the same as (13.61) after differentiation under 
the integral sign, resulting in the dependency of the cashflow upon the conditional 
payments EIly^u,t{Z) is solely through the rate of payment, meaning the derivative 
of t !->■ EIly^u,t{Z). As for (13.101) . this is simply the discounted expected payments 
after rearrangement of the order of integration. 
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For most policies of interest, the regularity conditions of Corollary 13.41 and Corol¬ 
lary 13.51 will be satisfied. Before proceeding to the next section, we give some 
examples of how to apply the corollaries in order to obtain expressions for cash¬ 
flows and liabilities for concrete types of spouse’s pensions. In Example 13.61 we 
derive a G82-type formula for a classical spouse’s pension, corresponding to the 
G82 reward no. 810. In Example 13.71 we calculate the liability of a terminating 
spouse’s pension, corresponding to the G82 reward no. 810, in a model where the 
health state of the spouse is modeled with longevity improvements. Finally, 
in Example 13.81 we consider a reward outside the scope of the G82 collection of 
rewards, where a lump sum is paid to the spouse at a particular age, and show 
how to use our results to derive the cashflow for this policy. 

Example 3.6. Recall that we in Example 13.11 identified the payment functions 
for a G82 reward no. 810 as given by the formula (13.3|) . We wish to calculate 
the liability of a reward no. 810 by applying Gorollarv 13.51 to this set of payment 
functions. To this end, we first need to specify a distribution for Z. Here, let 
be an inhomogeneous Markov chain on [u — j/, oo) with two states {a, d}, initial 
state a at time u and intensity r i—>■ qad{r + y — u) for transitioning from a to d and 
zero intensity for transitioning from d to a. This intensity corresponds to having 
Qadiy) be the intensity for dying at age y. Now, it is immediate from (|3.311 that 
^y,u,u{z) = 0 for all y,u>t). Furthermore, note that for y ^ 9, we have 

(3.11) = !(,>,) / P(Z“’^ = a)dr, 

J u 

so that 1 1 —>■ EIly u t{Z'^'y) is differentiable for t > u, and the derivative is 

= a) = exp qad{r + y - u) dr^ 

/ rv+t-u \ 

(3.12) =exp|^y qadir)drj, 

corresponding to the survival probability of the spouse from age y at time u to 
age y + t — u at time t. As a consequence, we obtain from Gorollarv 13.51 that the 
liability is given by 

(3.13) T= / h{u)g{u) / f{y\u) / e“’'‘exp[ / qad{r)dr\ dtdydu, 

Jo Jo Ju \Jy J 

corresponding to e.g. the formula for the liability found in Section 9 of the G82 
concession. o 

Example 3.7. In this example, we show how our generic model includes a G82- 
type model for a terminating spouse’s pension, corresponding to reward no. 815, 
including the application of longevity improvements when calculating the value of 
the terminating annuity for the spouse. A terminating spouse’s pension pays an 
annuity upon the death of the insured, conditional on the spouse having an age 
less than a given termination age c. As in Example 13.11 let E = {a, d}. Define, 
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for t > 0, 
(3.14) 


— ^(y^d)^(t>u) I l(j2r—o) l(y+r—?i<c) dr. 


The payment functions then correspond to the payments of a terminating 

spouse’s pension with termination age c. Next, assume given e.g. a nonnegative 
and continuous mapping q : K.+ x R+^R+, where q{t, y) represents the intensity 
for dying at time t when the age of the spouse is y at time t. Assume that Z*'’“ 
is an inhomogeneous Markov chain on [u — y, oo) with state space E, intensity 
r I— q{r, y-\-r—u) for transitioning from a to d and initial state a. This corresponds 
to that when the spouse has age y at time u, the intensity for transitioning from 
a to d at time r >u\s q{r, y + r — u). We then have Tly^u,u(z) = 0 for all y,u > 0, 
and for y ^ d, we have that t i—>■ u is differentiable for t > u, with 

derivative 


(3.15) 




= exp 


J qadir,y + r-u)dr^ , 


corresponding to the survival probability of the spouse from age y at time u to 
age y + t — u at time t. Corollary 13.41 then yields that the cashflow of the reward is 

(3.16) at = [ h{u)g{u) 

Jo 

corresponding to e.g. the formula for the liability found in Section 9 of the G82 
concession, with the exception of the intensity qad depending not only on age but 
also on time. o 


f{y\u)exp[ / qadir^y+ r-u)dr] dydu 


Example 3.8. In this example, we show how our generic model is capable of 
modeling non-G82 type rewards, in this case an insurance policy paying out a 
lump sum of one unit of monies to the spouse, when the spouse reaches age c, if 
the spouse at that time remains alive. If the age of the spouse at the death of 
the insured is c or older, the lump sum is paid out immediately. To define the 
payment functions, let E = {a, d}. Define, for t > 0, 

(3.17) — l(y/a)l(t>ii)l(y+t-u>c)l(zu+e-a=a)- 

These payment functions corresponds to receiving one unit of monies when achiev¬ 
ing the age c (corresponding to the indicator for y + t — u > c) foi any time t > u, 
or receiving one unit of monies at time u if the age c already is achieved at that 
time, with all payments contingent upon the spouse being alive at age c. Note 
that as c > 0 and z £ Du-y{E), the expression Zu+c-y is well-defined. As in 
Example 13.71 assume given a continuous mapping q : R_|_ x ]R_|_ —>■ ]R_|_, and assume 
that ^1'’“ is an inhomogeneous Markov chain on [u — y,oo) with state space E, 
intensity r i-A- q{r, y + y — u) for transitioning from a to d and having state a at 
time u. From Theorem 13.31 we know that the expected cumulative payments are 

pt poo 

(3.18) At= h{u)g{u) / f{y\u)Enu,y,t{Z'^^'^)dydu. 

Jo Jo 
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We wish to identify the corresponding cashflow. In contrast to the previous ex¬ 
amples, it does not generally hold in this case that ^y,u,u{z) = 0 for all > 0 , 
since lump sum payments can be made immediately upon the death of the insured. 
However, we do have that for y ^ 9 and t > u that 


(3.19) = ^y+t-u>c)PiZ::y,_y = a). 

Here, if c — y < 0, it holds with probability one that = a, since transitions 

from state c? to a are not possible and the state of Z'^'^ at time u is a. Therefore, 
we obtain for t > u that 


(3.20) EUy^uAZ'^^^) = My+t-u>c) = 1 for c < y, 

/ ru+c-y \ 

(3.21) i;n^,„,t(Z'“’^) = l(j^+t_„>c) exp ( y yad(r,y-|-r-M)drj for c > y. 

As a consequence, EIiy^u,t{Z'^’y ) is generally not differentiable as a function of t, 
since it may contain a jump at time c -I- u — y, corresponding to the time when the 
spouse reaches age c. Now define 

/ i<u+c-y \ 

(3.22) ^(u,y)=exp(y qad{r,y + r - u) drj dy, 


we then obtain for 0 < u < t that 


dt 


pOO pC 

f{y\u)Enu,y,t{Z^''^)dy = J f{y\u)dy + f{y\u)l(t>u+c-y)^{u,y)dy. 

f{y\u)^{u,y) dy 


As u + c — t < c, we obtain that 
d 


l0y{u+c—t) 


/(y|u)l(t>„+c_j,)C(M,y)dy = — 

(3.23) = l(u>t-c)f{c + u- t\u)^{u, c + u-t), 

where the derivative exists for t 7 ^ u -I- c. Furthermore, (13.201) and (I3.21|) yields 


(3.24) 


/ OO 

fiy\u)dy. 


Therefore, the Leibniz integration rule allows us to conclude that the cashflow 
exists and is given by 


at = 


dt 


K'a)g{u) / f{y\u)EXdu,y,t{Z'"'y)dydu 


poo 

+ h(t)g{t) f{y\u)Ent,y,t{Z^’y)dy 
Jo 

= J h(u)9{u)l(u>t-c)f{c + u-t\u)exp(^J qad{r,c + r - t) dr'^ dy du 

/ OO 

fiy\u) dy. 


Here, the first term corresponds to the expected payments per time given that 
the insured dies as married at time u < t. For payments to occur at time t, the 
spouse must have age c at time t, which corresponds to having age c + u — t at 
time u, which occurs with density /(c + u — t|u), and payments then occur when 
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the spouse survives to age c, which corresponds to the exponential. The second 
term results from the immediate payments made at time t when the spouse dies 
at t, and these payments occur precisely when the spouse has age c or more at 
time t, which corresponds to the integral of /(?/|u) over [c, oo). o 


4. An MPP model for marriage and spouse age with longevity 


Recall that the generic model outlined in Section [3] consists of variables T, X,Y 
and a set of processes for u,y > 0. In the expression for e.g. the expected 
cumulative cashflows, see (EH), the expressions involved are h(u), g{u), f{y\u) 
and EIlu,y,t{Z'^’y). In order to calculate cashflows and liabilities, it is necessary to 
model and calculate these four expressions. In this section, we develop a marked 
point process model allowing us to express g(u) and f{y\u) in terms of intensities 
for marriage, divorce and the death of the spouse. 

Note that specification of a model yielding a credible expression for h{u) is gener¬ 
ally not a problem, as this can be done by e.g. specifying a Markov chain model 
for the health state of the insured and letting T be the hitting time of the death 
state. Likewise, as we have seen in the examples of the previous section, in order 
to obtain explicit expressions for EIiu,y,t{Z'^’y), it suffices to specify e.g. a Markov 
chain model for the health state of the spouse conditionally on the time of death 
of the insured and the age of the insured at that time. 

In order to obtain a full specification of all expressions necessary to calculate 
expected cumulative cashflows et cetera, it therefore suffices to specify a model 
yielding expressions for g(u) and /(j/|m), which is the objective of this section. 
The motivations for considering this particular marked point process model of 
this section are twofold: We aim to obtain both a formalization of the model 
yielding the formulas of the G82 concession, as well as an extended model including 
longevity improvements for the spouse. 

In Subection l4.ll we construct the basic model framework, specifying a probability 
space with variables T, X and Y depending on an underlying marked point process. 
In Subsection l4.21 we derive expressions for g{u) and f{y\u) in this model. Finally, 
in Subsection 14.31 we show how to obtain the G82 expressions as a special case of 
our results. 


4.1. Model construction. In this subsection, we construct the marked point 
process based model for T, X and Y mentioned above. Our model must be such 
that we can express and calculate g(u) = P{X = Xm\T = u), the conditional 
probability of marriage given the time T of death of the insured, as well as the 
conditional density f{y\u) of F given T = u and X = Xm- As modeling of T is 
not our main interest, we will simply assume that a variable T with density h is 
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given. In order to obtain a joint model for T, X and Y, it remains to specify 
X and Y. To this end, we will develop a marked point process model for the 
combined marital state of the insured and the health state of the spouse, allowing 
for remarriage of the insured after the death of a spouse. We will assume that the 
distribution of the marital state of the insured at the time of death is the same as 
the distribution of the marital state at the time of death if the insured were not 
capable of dying. This removes the need for explicitly modeling the health state 
of the insured jointly with the marital state of the insured. 

We let denote the marital state of the insured at time t, with state space 
given by = {so,rni,si,m 2 ,S 2 , ■ ■ ■ ,}■ We think of the state Si as cor¬ 
responding to the i’th single (unmarried) state and think of as correspond¬ 
ing to the i’th married state. For convenience, we let s = {sqi Si, S 2 ,...} and 
m = {7711,7712,171.3,...}. Furthermore, we let denote the health state of the 
current spouse at time t, with state space E'^ = {a,d,d}, where, similarly to the 
previous section, 3 is a coffin state, in this case the state held by the process Z^ 
when the insured is unmarried, a corresponds to the spouse being alive and d cor¬ 
responds to the spouse being dead. Finally, we let denote the age of the spouse 
at time t, with state space K+ U {i9}. We assume that {Z^, Z‘’,U‘’) is independent 
of T. 

Our next task is to construct and specify the joint distribution of the processes 
{Z^, Z^, U^) in the model. To this end, first define stopping times 

(4.1) Tt = inf{t > 0 I = 1 } 

for any i G E^^. Assume given an uncountable sequence of variables (Y^’^)t>o, all 
independent. We interpret T'”’* as the hypothetical initial age of a spouse married 
at age t. We assume that the variables Y^ are such that Y^ has density </ 7 (-|t). 
We then define 

(4.2) = +(t-T,0J 

(4.3) Uf = Uf'‘' when {Zl^,Z^) G {( 77 ii,,a), (sj,,a)} and = d otherwise. 

Note that in (14.21) . the process increases indefinitely after Tj)) . Consistently 
with this, is the hypothetical age of the u’th spouse given that the spouse 
cannot die. We then let {Z^,Z‘’) be a marked point process with state space 
E^ X E^, initial state (s,i9) and with intensity A : R+ x {E^ x E‘’) x {E'^ x E‘’) 
given by 

(4.4) A(t, (sj._i,i), {m„,aj) = 7 (f) for all i G E'^ , 

(4.5) \{t,{my,a),{sv,a)) = a{t), 

(4.6) \{t, dj) = g'^(t, I/f), 

and all other intensities zero. For details on marked point processes and intensities, 
see [7]. Here, 7 : K+ M- 1R+ denotes the marriage intensity, a : 1R+ 1 —>■ K+ denotes 
the divorce intensity and : 1 R+ x 1 R+ 1 —>■ ]R_|_ denotes the death intensity for 
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the spouse, with q^{t,y) denoting the intensity at time t when the spouse has 
age y at time t. Note in particular that the states (so, a) and (so,d) never occur, 
corresponding to the impossibility of having a spouse in the initial single state Sq. 
See Figure [47T] for an illustration of the possible transitions of the process. 



Figure 4.1. Transition graph of the combined marked point pro¬ 
cess for the marital state of the insured and the health state 
Z^ of the spouse. 


We have now made sufficient assumptions on the processes Z^ ,Z‘’ and C/^ to 
uniquely determine their joint distribution. We then simply define 


(4.7) 

(4.8) 


X — I when Z^ G s 

\ Xm when Z!^ G m 

Y = U^. 


This yields a complete model for T, X and Y, and we may then express g(u) and 
fiylu) in terms of the underlying marked point process. To see how this can be 
done, note that by the independence of T and (Z^, Z^,U^), we have 

(4.9) g(u) = P(X = Xm\T = u) = P(Z!^ G m\T = u) = P{Z!^ G m). 

Furthermore, for measurable A C R_|_, we have 



u) dy = P{Y gA\T = u,X = Xm) = P(U^ £ A\T 


u, Z!^ G m) 


(4.10) = P(Ui G A\T = u,Z!^Gm)= P{U< G A\Z^ G m). 


Since A was arbitrary, this yields that f{-\u) is the density of given Z!^ G m, 
assuming that this density exists. In conclusion, we now have that 

(4.11) g(u) = PiZ): G m) 

(4.12) /(» = -^P{Ui G ■\Z>^ G m), 

where denotes the Lebesgue measure and the derivative refers to the Radon- 
Nikodym derivative. Our next objective is to obtain expressions for (|4.11|) and 

(SH. 
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4.2. Expressions for the marriage probability and spouse’s age density. 

In this subsection, we state expressions for (14.1111 and (|4.12l) in terms of the inten¬ 
sities for marriage, divorce and spouse’s death. To do so, we first introduce the 
auxiliary expressions 


(4.13) u^t) = PiZt = s,) for u = 0,1,..., 

(4.14) 5i^(-|0 = e •, for u = 1, 2,..., 


and further define 

(4.15) = exp 

(4.16) = exp 

(4.17) = exp 


J 7 (u) dv^ , 


a{v) du 


/ q^iv,y + v 

Jo 



The interpretation of the latter expression is that i^y is the survival probability 
for a spouse starting at time zero and until time t, with the assumption that the 
spouse, if surviving, will be y years at age t. In Theorem 14.41 and Theorem 14.51 
we state results on how to express (j4.11[l and (14.1211 in terms of these auxiliary 
expressions. Before stating the two theorems, we first state three lemmas, yielding 
expressions for (I4.13|l and (j4.14ll . 


Lemma 4.1. It holds that uo(t) = €[ for t > 0. 


The following lemma shows that the density defined in (I4.14|l in fact exists, and 
yields an expression for the density. In the following, we let denote the Borel 
(T-algebra on K+. 


Lemma 4.2. Let u > 1 and define a measure Qt,v '■ B+ — >■ [0,1] by 
(4.18) QtAA) = P{U^ e 41, Zf = m,). 


Then has a density with respect to the Lebesgue measure, and the 

density is given by 


(4.19) 


9v{y\t)=l u,j-i{vA{v)Ay -t\v)- 




^v,y+v—t 


du 


Lemma 4.3. For v >1, it holds that 

nCO pt 

(4.20) u 


nao nt £7 

it)= j / g{y\v){cr{v)+q^{v,y))^dvdy. 
Jo Jo 


Note that combining Lemma 14.11 Lemma 14.21 and Lemma 14.31 we obtain expres¬ 
sions for Uv{x) and gu(ri\x) in terms of coupled recursion equations. The inter¬ 
pretation of the above lemmas is as follows. As regards Lemma 14.11 this lemma 
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simply states that uo{t) is the survival probability for the distribution with haz¬ 
ard 7 , corresponding to that the only way of leaving the first single state so is to 
become married with intensity 7 . The formula (|4.19p expresses that the density of 
being married for the iz’th time at time t with a spouse of age y at time t can be 
obtained by, for some v < t, being single for the iz — I’th time at time v, becoming 
married with intensity 7 (u), to a spouse which at time t has age y, corresponding 
to having age y + v — t at time u, and finally, not being divorced in the time 
interval [u,t] with intensity cr and not having the spouse die in the time interval 
[u,t] with intensity r i-)- q‘’{r,y + r — t). Similarly, the formula (14.201) expresses 
that the probability of being single for the z^’th time at time t, for v > 1, can be 
obtained by conditioning on being married for the i/’the time at time v < t with 
a spouse whose age at time v is y, occurring with density 5 i/(i/|u), and afterwards 
either being divorced with intensity <j(v) or having the spouse die with intensity 
q‘’{v,y), and finally not remarrying, where marriage has intensity 7 . 

With these lemmas at hand, we are now ready to state our main results yielding 
expressions for g{t) and f{y\t). 

Theorem 4.4. With g{t) = P{Z^ G m) as in {4.11^ , it holds that 

00 pQQ 

(4-21) 5W=X!/ 9v{y\t)<^V- 

Theorem 4.5. Define Q* : S+ i->- [0,1] by Qt{A) = P{U^ G A\Zf = m). It 
then holds that Qt has a density with respect to the Lehesgue measure. With f(-\t) 
denoting the density as in the density is 

00 

(4-22) f{y\t) = -77^Y1 

git) ^ 


Theorem 14.41 and Theorem 14.51 expresses g{t) and f{y\t) in terms of the auxiliary 
density giy{y\t), which in turn is obtained from the fundamental model parameters 
and Uyff). Essentially, (14.211) states that the probability of being married can 
be obtained by integration of the density corresponding to being married with 
a spouse of a particular age, and (14.221) states that the density of the spouse’s 
age given marriage can be obtained by considering the density of the measure 
corresponding to being married with a spouse of a particular age and normalizing 
with the probability of being married in order to obtain the conditional density. 


4.3. Comparison to the G82 model. In this subsection, we show how to re¬ 
claim the expressions for the marriage probability and spouse age density from 
the G82 model, see Section 8 of the G82 concession. To this end, consider the 
model setup of Subsection l4.ll and assume that the insured at time t has age zero. 
We further assume no longevity improvements for the spouse, meaning that we 
let q^(t,y) = q''{y), with minor abuse of notation, corresponding to q‘’{y) denot¬ 
ing the intensity for the spouse’s death at age y. Furthermore, we assume that 
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a{t) = 'y(t) = 0 for t < a, corresponding to a lower age bound for marriage. In 
order to obtain consistency with the G82 formulas, we let x denote time in the 
following, or equivalently, the age of the insured, and we let rj denote a generic 
value of the age of the spouse. Furthermore, we extend g^(r]\x) and f(ji\x) from 
K+ to K by letting gvivlx) = fig\x) = 0 for 77 < 0 . 


With this notation, Theorem 14.41 and Theorem 14.51 yields 

00 pOo 


(4.23) 

(4.24) 


pOO 

9{x)=^j gv{'n\x)(l'q, 

,,_1 J —CO 


' —CO 

00 


9ix) ^ 


Now define 

(4.25) = exp J q‘^{y) dy^ . 
Recalling (14.171) . we then have 

= exp q‘'{v,ri + V — x) dv^ 

(4.26) = exp q'^ {f] + V — x) dv^ = 


f 

^x,r\ 




^TJ+^-X 


Therefore, Lemma [4.21 and our assumption that 7 (^) = 0 for ^ < a yields 


9Av\x)= [ i - x\0^ J"'’'' dg 

(4-27) =f + V - x\0^ 


Furthermore, Lemma 14.31 yields 


u^{x) 


(4.28) 


/ OO pX py 

/ 9{y\0{<x{0+ 

-00 Jo 

/ OO pX ^7 

/ 9ii + g- a:|C)(cr(C) + qH^ + V-x))-^d^ dg. 

-OO Jo 


Here, the change of variable from y to ^ + g — x corresponds to a change of variable 
from the spouse age y when the insured has age ^ to the spouse age g when the 
insured has age x. Finally, the assumption that 7(0 = 0 for ^ < a immediately 
yields that uo(a;) = 1 for a: < a and uo{x) = tll^l for x > a. All the above 
formulas now match exactly the formulas from Section 8 of the G82 concession, 
showing how to obtain these formulas as a special case of our generic model in 
general and our marked point process model in particular. 
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5. Discussion 


In this paper, we have introduced a generic model for spouse’s pensions, and have 
derived expressions for expected cumulative payments, cashflows and liabilities. 
By example, we have shown that these results can be used for the calculation of 
the cashflows and liabilities for both classical G82-type spouse’s pensions as well as 
for more advanced policies such as the policy considered in e g. Example l3.81 Fur¬ 
thermore, we have developed an explicit joint marked point process model for the 
marital state of the insured and the health state of the spouse, including longevity 
improvements for the spouse, and have shown how to derive the expressions rele¬ 
vant for cashflows et cetera using this model. Finally, we have shown that in the 
absence of longevity effects for the spouse’s health state, the expressions for the 
marital probability and the density of the spouse’s age from the G82 concession 
can be reclaimed from the results obtained here. 

The immediate benefit of the generic model developed in Section [3] is the result 
that common expressions for cashflows and liabilities hold independently of the 
particular model for marital probabilities and spouse’s ages. This allows for a 
simple framework for deriving liabilities for various types of spouse’s pensions 
without consideration of the particular model for marital behaviour to be applied. 

The results also makes explicit that one main issue for the calculation of liabilities 
for spouse’s pensions is the estimation of the probability of marriage g{t) and the 
density of the spouse’s age f{y\t). One opportunity for further work is to consider 
methods for efficient estimation of these functions, e.g. through direct estimation 
from the marital states of a general population, or through estimation of intensities 
in a more specific model such as the marked point process model considered in 
Section m 

The results of Section 0] opens up for the possibility of obtaining more accurate 
estimates of liabilities by taking into account the longevity improvements of the 
spouse when calculating the probability of marriage for the insured. We take a 
moment to reflect upon the complications of numerical computations resulting 
from this extension. In the classical G82 model, one component of the numerical 
effort required for computations is the evaluation of g{x) and /{rilx), see Subsection 
14.31 for X = 0,..., 125 and rj = 0,..., 125, with 125 conventionally being taken 
in the G82 concession as the maximal age of the insured and the spouse. In 
Subsection 14.21 these expressions are parameterized in terms of time, i. e. g(t) 
and f{y\t), and depend on given intensities 7 and cr for marriage and divorce, 
respectively. For the liabilities corresponding to a single insured, therefore, the 
computational effort required is no greater than in the case without longevity. 
For an entire population, however, some computational overhead occurs: The 
intensities 7 and cr correspond to marriage and divorce intensities as a function 
of time. Given a model where these intensities ultimately depend on the age of 
the insured, we would have e.g. j{t) = 7 a(a:o + 1) and a{t) = aa(xo + t), where 
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7 a and Ua denote the intensities for marriage and divorce as functions of the age 
of the insured, and xq denotes the age at time zero of the insured. As a result, 
the marital expressions g(t) and /(?/|t) would in fact depend on the initial age 
of the insured. For the calculation of the total liabilities for a pension fund, it 
would therefore generally be necessary to calculate a grid of values g{x,xo) and 
/(r 7 |a;,a;o) for a;o = 0,..., 125, x > xq and y = 0,..., 125. Other than this, 
no particular increase in computational complexity would be incurred from the 
inclusion of longevity improvements. 

As regards opportunities for further work, the most pressing necessity for the 
accurate estimation of liabilities is the estimation of g{t) and f{y\t). The issue 
of obtaining standard methodology for this, either by direct estimation or by use 
of more complex models such as the marked point process model, is not very 
well developed. Furthermore, as the calculation of g{t) and f(jj\t) in the marked 
point process model is relatively computationally intensive, it is of interest to 
develop models for the marital behavior of the insured and the spouse which is 
both amenable to estimation of parameters and in which simple expressions for 
g{t) and f{y\t) can be obtained. 


Appendix A. Proofs 


This appendix contains proofs of the results of the main part of the paper. 

A.l. Proofs of results in Section [3l In this subsection, we prove the results of 
Section[3]on expected cumulative payments, cashflows and liabilities in the generic 
model for spouse’s pensions. 

Proof of Theorem 13.31 Fix t > 0 and let ikt : K+ x {xs,Xm} x FV —>■ M be 
defined by 



(A.l) 


Recalling (13.2L it then holds that Bt = '^t{T,X,C). With R denoting the dis¬ 
tribution of {T,X) and (Qu,x) denoting the conditional distribution of C given 
{T,Xm), we then obtain that 


EBt 


E^t{T,X,C)= [ 

J(i 


{R+U{0})x{a::s,aim} - - 



J 



(A.2) = [ f A(t)-A(u-)dQa.,(A)dR(u,x) 


J[0,t]x{£Cm} JFV 
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Now note that 


i?([0,t] X {x^}) = P{T <t,X = Xm) 

(A.3) = f h{u)P{X = Xm\T = u) du = 

Jo 

so that with tt : IR.+ x {xm} R+ defining the projection mapping onto the 
first coordinate, it holds that the pushforward measure of the restriction of R 
to IR.+ X {xm} under tt has density u i—>• h{u)g(u) with respect to the Lebesgue 
measure. Inserting this into the above yields 

(A.4) EBt= [ h{u)g{u) [ F{t) - F{u-) dQu,xmiE) du. 

Jo ^FV 

Furthermore, by our assumptions, the conditional distribution of the variable Y 
given T = u and X = Xm has density f{-\u). Therefore, the conditional distribu¬ 
tion of (T, Y) given T = u and X = Xm is the tensor product of the Dirac measure 
in u and the measure with Lebesgue density f{-\u). With Qu,x,y denoting the 
conditional distribution of C given T = u, X = x and Y = y, we therefore obtain 
for measurable A C FV that 


h{u)g{u) du, 


(A.5) 



Qu,x,n,y{^)fiy\'^) ^y- 


Next, recalling that Ct = ), and further recalling that for all 

u, y > 0 is assumed to be independent of (T, X, Y), we obtain that the conditional 
distribution Qu,xm,y is equal to the distribution of the process t i-A- Uu,y,t{Z‘^’'^), 
we obtain 


[ F{t)-F{u-)dQu,xAF) 

Jfv 

= [ fiy\u) [ F{t) - F{u-)dQu,xm,yiF)<Jy 

Jo Jfv 


(A.6) 




Collecting our conclusions and recalling that Ilu^y^u-iz) = 0 for all z G Du{E), we 
finally obtain 

nt poo 

(A.7) EBt= h{u)giu) /(y|u)An„.y,t(Z“'^) dy du, 

Jo Jo 


as required. 


□ 


Proof of Corollary 13.41 By Theorem 13.31 we have 


(A.8) 


EB, = 


Hu)g{u) / /(y|u)FTIj^_„^t(Z“’^) dydu. 
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Applying the Leibniz integral rule, we then obtain 

(A.9) at = -^EBt = J ■^h{u)g{u) J f{y\u)EIly^u,t{Z^’'^) dy du 

pOO 

(A.IO) +h{t)g{t) / f{y\t)EUy,tAZ^’y)dy. 

Jo 

By our assumptions, the latter term is zero. Next, as y f{y,u) is a probabil¬ 
ity density, any bounded interval is integrable with respect to the measure with 
Lebesgue density f{y\u). As we have assumed that EIly^u,t{Z^’y) is bounded as a 
function of t on compact intervals, this yields that differentiation under the inner 
integral in (IA.9I) is allowed, and we obtain the result stated in the corollary. □ 


Proof of Corollary 13.51 By approximation with Riemann sums and a uniform 
integrability argument, we have 


(A.11) 


L = E 


e-^UBt = 


f 


e’'* dEBt. 


Applying Corollary 13.41 we then obtain 


L = 

(A.12) 
as required. 


/ ^ j j f(y\'^')(^'^ZIly^u,t{Z)^ dydudt 


h[u)g{u) / f{y\u) / e-*-* 


—EBy^u,t{Z) ) dtdydu. 


□ 


A.2. Proofs of results in Section In this subsection, we prove the results 
of Section 0] on expressions for the marriage probability and the density of the 
spouse’s age. 

Proof of Lemma 14.11 This follows as 

(A. 13) uo (0 = P{Zt = So) = P{Zl^ = so,Z^ = d) = exp J j{v) du^ , 

since the intensity for leaving state (sq; d) is 7 . □ 

For the following lemmas, we require some results on compensators and intensi¬ 
ties. For general results on compensators, see Chapter V of [3]. For results on 
compensators and intensities in the particular context of marked point processes, 
see Chapter 3 of [7]. 

Lemma A.l. It holds that has a density : R+ — >■ R+ with respect to the 
Lebesgue measure, and the density is given by hra,,{t) = PiZ^ = s,^-i)y{t). 
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Proof of Lemma lA.il Define A™*' by 
(^•14) A™" = ^ 

0<v<t 


It then holds that the compensator A^'' of A™*' is 

(A.15) A^’' = f -f{v)l(^z>;=s,_^)Av- 

Jo 

Since only a single jump to is possible, this yields 


<t) = ENr = =E f 7(^^)l(zr dv 

Jo 

(A.16) = [ P{Z^ = Sj,-i)^{v)dv, 

Jo 


proving the lemma. 


□ 


Proof of Lemma 14.21 Applying Lemma lA.ll and (Z^ = m^) C < t), we 
obtain 


P{U^ e A, = m,) = P{U^ e A, Zt = < t) 

P{U^ G A, Z^ = = v)P{Z^ = s^.Oliv) Av 

(A.17) =[ u^-i{v)j{v)P{U^ € A,Z^ = m^\TI^^ =v)dv. 

Jo 

The intensity of Z^ leaving state TOj, at time u is a{u) + q‘’{u, U^). Therefore, we 
obtain 



P{U^ € A, Z^ = m,\TI},^=v) 

= P{Y< + {t-v)€ A, Zt = = v) 

(A.18) = f ip{y + V - t\v)^ dy. 

From this, we obtain that P{Uf G •, Z^ = = v) has a density given by 

(A.19) y ^ ip{y + V - t\v)^ . 

Combining (|A.17p and (IA.19p . we find that Qt^^, has a density, and the density is 

1 pt 

9i^iy\t) = “T / < y, = m^\T^ = v) dr 

dy Jo 

= [ u^-i{v)-f{v)-^P{U^ <y,Z^ = m„\TI^^ =v)dv 
Jo dy 

(A.20) = [ u^-i{v)-f{v)qy(y + v-t\v)^ dr, 

Jo V 
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as required. □ 

Lemma A.2. Let u > 1. It then holds that Tjl has a density hg^ : M+ —>■ R+ with 
respect to the Lebesgue measure, and the density is given by 

(A.21) hs,,{t) = g„{y\t)(a{t) + q‘^{t,y))dy 

Jo 


Proof of Lemma IA.2i Define a process by putting 

(^•22) = X! 

0<y<x 


It then holds that the compensator of N‘^'' is 


(A.23) = / l(z?=m,)(o-(w)+g‘’(w, C/5))di;, 

Jo 

which yields 

(A.24) P(T't <t)= f + qHv, U^)) 

Jo 


so that has density given by hs,,{t) = It remains 

to show that the density can be written on the form (IA.2II) . In order to obtain 
this, let Qt^i, be the measure defined in (j4.18|) . and recall that by Lemma Qt^i, 
has a density given by (14.191) . Note that for all A S P, it holds that 


(A.25) 



l(^=m„)lA(2/) d{z^, U^){P){z, y) 


1 a ( 2 /) dQt,,y{y), 


with {Z^ ,Uf){P) denoting the distribution of {Zf,Uf). As a consequence, the 
same relationship holds with Ia('w) exchanged with arbitrary measurable and 
bounded / : R —>• K. Applying this, we obtain 


(A.26) 


hsAt) = El(^z^=m,,){^{t) + 

= f l(z=m„)(o-(t) + g‘^(t,j/))d(Zf,[/t‘^)(P)(z,2/) 

JE^^xR 

= f cr(t) + qHt,y)dQt,i.i.y) = f gu{y\t){cr{t) + qHt,y))dy, 
Jr Jo 


as was to be shown. 


□ 
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Proof of Lemma I4.3L Noting that (Zf = s,y) C (Tj^ < t) and applying Lemma 
IA.21 we obtain 


.(t) = P(zr = s,) = P(Zf = < t) 

pt poo 

= PiZt=Su\T^^=v) g{y\v){a{v)+q^{v,y))dydv 
Jo Jo 

poo pt 

g{y\v){a{v) + q^{y,y))P{Z’^ = = v)dvdy 


/o ^0 

nOO pt 


(A.27) 

as required. 


/o ./o 


fi 

g{y\v){a{v) + q^{v,y))-P dv dy, 

Zv 


□ 


Proof of Theorem 14.41 Applying the notation and results of Lemma 14.21 we 
obtain 

OO OO pOO 

(A.28) g(t) = P{Z^ ern) = ^ = X! / 9u(y\^) 

u=l 1 /= 1"'0 


as required. 


□ 


Proof of Theorem 14.51 Recalling Lemma 14.21 it holds for A G B+ that 
P{U^ G A, Z>^ = m) 1 


QM) = 


(A.29) 


= m) g[t) 

OO n ^ 


J2PiU} € A, Z^ = m) 


1^ = 1 


git) 


'^g„iy\t) dy, 




implying the result. 


□ 
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